Resonant plasmon scattering by discrete breathers in Josephson junction ladders 
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We study the resonant scattering of plasmons (hnear waves) by discrete breather excitations 
in Josephson junction ladders. We predict the existence of Fano resonances, and find them by 
computing the resonant vanishing of the transmission coefflcient. We propose an experimental setup 
of detecting these resonances, and conduct numerical simulations which demonstrate the possibility 
to observe Fano resonances in the plasmon scattering by discrete breathers in Josephson junction 
ladders. 
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INTRODUCTION 

Discrete nonlinear Hamiltonian systems generically al- 
low for spatially localized and time periodic states, dis- 
crete breathers (DBs), which exist thanks to the inter- 
play between nonlinearity and discreteness Q . DBs have 
been detected and studied experimentally in interacting 
Josephson junction systems coupled nonlinear opti- 
cal waveguides Q , lattice vibrations in crystals , anti- 
ferromagnetic structures 0, micromechanical cantilever 
arrays [a, Bose-Einstein condensates loaded on optical 
lattices pi, layered high- Tc superconductors j^. DBs are 
predicted also to exist in the dynamics of dusty plasma 
crystals 0. 

In the transmission problem of small amplitude waves 
at frequency through a DB in the particular case of 
a one-dimensional lattice, the DB acts as a time peri- 
odic scattering potential with frequency Q. Generally, 
there is an infinite number of scattering paths in such a 
case, which may lead to a variety of interference phenom- 
ena (see Fig. QJ. Indeed, perfect reflection was observed 
for particular wave numbers [lol [ill flp . IT^ . A detailed 
analysis of this phenomenon shows [14| that it is a Fano 
resonance 15] which is based on the phenomenon of de- 
structive interference. 

Systems of coupled Josephson junctions with a ladder 
geometry allow for the experimental excitation and de- 
tection of DBs A variety of experimental methods 
allows for a systematic study of DB properties, which 
have been successfully combined with theoretical calcu- 
lations and predictions [3, [13, El El iSa, 21]. JJLs 
are quasi-one-dimensional, and allow for the excitation 
of travelling linear oscillatory waves (usually referred to 
as plasma waves or plasmons), which can be scattered 
by DBs. That makes JJLs suitable for the observation 
of Fano resonances. In addition, Josephson junction sys- 
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FIG. 1: Schematic representation of dynamically generated 
paths in the frequency domain by a time periodic scattering 
potential with a period 2-k/Q.. In general, there is an infi- 
nite number of paths which may lead to different interference 
phenomena. 



terns are dissipative systems including both damping and 
the presence of an external homogeneous dc bias. 

The aim of this paper is to investigate the possibility 
for the observation of Fano resonances in plasma wave 
scattering by DBs in Josephson junction ladders (JJLs). 
We will also discuss the experimental setup for the possi- 
ble observation of these resonances. We make use of the 
underdamped regime of dissipation and exploit previous 
results on Fano resonances which have been obtained for 
nondissipative systems. 
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MODEL 



JJLs are formed by an array of small Josephson junc- 
tions that are arranged along the spars and rungs of a 
ladder, as shown in Fig. [21 Junctions are denoted by 
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FIG. 2: Schematic view of the Josephson junction ladder. 
Each cross represents a junction. Arrows denote the dc bias 
flow. 

crosses. Each junction consists of two small weakly cou- 
pled superconducting islands. The dynamical state of a 
junction is described by the phase difference (j) (Joseph- 
son phase) of the superconducting order parameters of 
the two islands. When the phase difference does not vary 
in time 4> = const, the junction is in the superconduct- 
ing state. Otherwise, the junction is in a resistive state 
with a nonzero voltage drop V a (j) (the dot represents 
differentiation with respect to time). 

By using the resistively shunted junction (RSJ) model 
and Kirchhoff's laws, we obtain a set of equations for 
Josephson junction ladders (JJLs) in the following form 

C + aC + sinC = 7+-^(AC + V</.ti-V^ti) 

PL 

^i^j-f «<^:j-Ksin</.:j = -±-{vci,i+4t-4,'^^) (1) 

^f, H- + sin = ^(VC + 4>i - ~4>i) , 

VPL 

where the notations A/„ = fn-i — '2fn + fn+i and 
V/„ = fn+i — fn are used. The Josephson phases of 
vertical, upper and lower horizontal junctions in the nth 
cell are denoted by ^J^, (j)'^ and (p!^ respectively, a is an 
effective damping, 7 is the external dc bias across vertical 
junctions. The discreteness parameter Pl = 27r/^L/<I>o 
characterizes the ratio of geometrical cell inductance L 
and the Josephson inductance of vertical junctions, and 
ry = /Ic is the anisotropy parameter of the Ladder. 
Here, lY {I^) is the vertical (horizontal) junction critical 
current (see Ref. for details). 

SMALL AMPLITUDE EXCITATIONS OF THE 
SUPERCONDUCTING GROUND STATE 

First we discuss the spectrum of small amplitude ex- 
citations (plasmons) around the superconducting ground 



state c^*" = arcsin7 and 0*'' = <j)^ = by linearizing 
the system ^ 

<\>l = 0r + , € = + A , 0' = + , (2) 

where ip"^, ip'^^ and (f'^ describe the small amplitude ex- 
citations. The presence of dissipation leads to the decay 
of excitations in time. In experiments, the dissipation a 
can be rather weak (of the order of 0.01 or even less) and 
usually it is ignored when discussing the small amplitude 
excitation spectrum ,161] . Here we provide the correct so- 
lution in the presence of damping. We can exclude the 
dissipative term by the following time-dependent trans- 
formation 

y) = e"^*V'- (3) 

According to Q the obtained plane waves are charac- 
terized by an exponential decay in time with a char- 
acteristic decay time t = 2 /a. By solving the equa- 
tions for V'n ~ Q^ii^'^-'^qt) obtain three plasmon 
bands LUi{q) , uj2{ci) , ^3{q) (see Appendix A). One of them 
Luf{q) = l — ^ is dispersionless (i.e. the frequency does 
not depend on q). It is characterized by all vertical junc- 
tions being at rest and in-phase excitations of horizontal 
junctions in each cell. The two other plasmon bands are 
characterized by a;|(g) < Wi(q) and uj'^{q) > uji{q)- Note 
here, that the band 102(1) possesses a rather weak disper- 
sion and becomes dispersionless for 7 = 0. In that limit 
it also coincides with the first band: ^2(9) — ^i('z) for 
7 = 0. Therefore, in the Hamiltonian limit a = and 
7 = there is only one plasmon band wl (q) with nonzero 
dispersion. 

DISCRETE BREATHERS 

JJLs support dynamic localized states - discrete 
breathers. A breather is characterized by a few junc- 
tions being in the resistive state (</>) while the others 
reside in the superconducting state {(p) = 0. The fre- 
quency of a DB is proportional to the average voltage 
drop across the resistive junctions fib oc (0) and gener- 
ally depends on the parameters of the system. Different 
types of DBs have been observed experimentally and nu- 
merically |jl6|. In the following we focus on DB solutions, 
which are schematically represented in Fig. |31 

By tuning the external dc bias 7, a DB state gener- 
ally follows adiabatically and change its characteristics 
including voltage drop (frequency), spatial extent etc. 
We coin such a state a nonresonant DB. This is correct 
unless the DB solution starts to resonate with the plas- 
mon modes of the unexcited part of the ladder. Such 
resonances may either lead to a loss of the DB solution, 
or to the appearance of resonant DBs [Hill 113. These 
states are characterized by a strong coupling of the DB 
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FIG. 3: Schematic representations of A- and B- type breather 
resistive cores in JJLs with one vertical resistive junction. 
Black spots mark the junctions being in the resistive state 
while all others are in the superconducting one. 

and the extended plasmon modes of the ladder, and a cor- 
responding voltage locking. Consequently resonant DBs 
are characterized by their voltage drop (frequency) be- 
ing nearly independent of the applied external dc bias. 
Note that resonant DB states do not persist in the dis- 
sipationless Hamiltonian limit easily, since the absence 
of dissipation causes a radiation of the DB energy into 
the rest of the ladder due to the assumed resonance with 
the extended ladder modes. In the dissipative case, ex- 
tended modes are damped, and energy is fed into the 
system through the breather core, providing with an in- 
tricate balance. However DBs with nondecaying tails in 
Hamiltonian systems might well correspond to resonant 
DBs in dissipative systems. 

DISSIPATIONLESS LIMIT AND 
CORRESPONDENCE IDEA 

Let us study the Hamiltonian limit, by neglecting both 
the dissipation a and the dc bias 7. In that case we can 
later apply a well developed numerical scheme for wave 
scattering. While the wave scattering can be also ana- 
lyzed for the dissipative case, we use the advantages of 
DB properties in Hamiltonian systems. In this limit non- 
resonant DBs become one-parameter families of periodic 
orbits in phase space. Their general form is 

$n{t) = Kn^bt + gn{t) , (4) 

where the notation /„ = {/^, /^} for the function of 
cell n is used. The jth component of Kn is an integer 
winding number. It is nonzero for a few junctions in the 
DB's core and zero for all others. E.g. for the breathers 
from Fig. 13 we have k = 1 for all resistive junctions of the 
A-breather and zero otherwise, while for the B-breather 
K — 1 for the resistive horizontal junctions and k = 2 for 
the vertical resistive junction. gn{t) are periodic func- 
tions of time gn{t + ^) = gn{t) with exponential decay 
of amplitudes along the ladder max|g? , 0. 



The frequency fib is a parameter which changes along 
the family of solutions. 

In the presence of weak dissipation a and dc bias 7 a 
given solution from a nonresonant DB family is selected 
to become an attractor and only slightly changed in its 
temporal evolution. These changes are small if a <C 1 
which we assume. By tuning the external dc bias 7 we 
simply scan dissipative nonresonant DBs which are close 
to their nonresonant DBs of the Hamiltonian limit. For 
this purpose one can use the relation between f2f,, a and 
7, which is given in Ref. [l^ . This correspondence idea 
will be exploited in what follows below. 

SCATTERING OF WAVES BY DISCRETE 
BREATHERS: DISSIPATIONLESS CASE 

DBs are dynamic localized excitations. For propagat- 
ing waves DBs act as time-periodic scattering potentials. 
Based on results for Hamiltonian lattices JOj, JJ^, , 
we expect that such interesting phenomena as resonant 
transmission and resonant reflection could be observed 
in JJLs too. Here we focus on the resonant reflection, 
or Fano resonance. We start with the Hamiltonian limit, 
where we can rely on recent results on Fano resonances 
induced by time periodic scattering potentials 14]. We 
predict the position of the resonance and then conflrm the 
prediction using numerical simulations. The extrapola- 
tion of the obtained resonances to the non-Hamiltonian 
case in the presence of non zero damping and bias will 
be shown to be successfull. 

In order to study the scattering problem in the Hamil- 
tonian limit a — and 7 = 0, we linearize the set of equa- 
tions m around a DB Q by substituting ^„ = 0„ -f 0^ 
with small amplitude excitations ipn 

PL 

(p^ + COsi^t)^^^ = + - ^n) (5) 

This is a set of linear equations with time-periodic coef- 
ficients with the period Tb = jf'- 

The set of equations (O describes scattering of waves 
by a time-periodic scattering potential. We may esti- 
mate the time-averaged scattering potential by replacing 
cos{(j)) = for a resistive junction and cos(0) = 1 for 
a superconducting junction. The interaction part (right 
hand side of Eq.jSl) is not changed. Due to the finite 
maximum strength of the scattering potential, the ex- 
pected nonresonant scattering impact can be expected to 
be weak in general, especially if the plasmon band width 
(and thus the corresponding group velocity of waves) is 
large enough. That promises ideal grounds for observing 
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resonant reflection on the background of nearly perfect 
transmission. 

The wave scattering by DBs is studied by using proper 
boundary conditions, which implies a scattering setup - 
at the left hand side there are incoming and reflected 
waves and at the right hand side there are only transmit- 
ted ones 

^ J /e-Ut^a (<?)«-«") ^^e-i(c^3(9)t+<?n)^ " < /R^ 

with the frequency ^3(5) H21|l and relation (|22(l between 
the components of the vectors /, R and T. 
The general solution of lO can be written as 



-i{uJ3{q) + kQt)t 



(7) 



Each term in this sum represents a channel - the path 
way for the waves. There is an infinite number of chan- 
nels. The fcth channel is characterized by its frequency 
u)g = u}3{q) + kVLb. When the frequency w^' belongs to 
the spectrum (|21|l the fcth channel becomes open, other- 
wise it is closed. It leads to the following property for the 
amplitudes B^n 



B 



k,\r 



= 0, for closed channels . 
7^ 0, for open channels 



(8) 



In other words, plane waves can freely propagate only 
inside the open channels. 

The zeroth channel fc = is always open. In principle 
there is a possibility to have an additional open channel 
uj™ with some nonzero odd integer m =/= implying 

uj3{q) + ujs,iq') = -mflb . (9) 
This multi-channel scattering takes place when 



\^3{0)\ < < k3(7r)| . 



(10) 



It corresponds to the case of parametric resonance, which 
may lead to the instability of a DB. Moreover, it was 
shown that multi-channel scattering is an inelastic pro- 
cess III the following we restrict our attention only 
to the case of one-channel scattering. We can avoid the 
multi-channel scattering by choosing frequencies of DBs 
such that the inequality 110|l is not satisfied. 

In the case of one-channel scattering only the zeroth 
channel fc = is open and all other channels fc ^ are 
closed. From this point of view the scattering potential 
which is generated by the DB, can be separated onto two 
parts: a time averaged dc part (the open channel) and 
an ac part (the closed channels). The scattering by the 
dc part only is well understood. We are interested in the 
role of the set of closed channels which become active in- 
side the breather core and may provide with interference 
phenomena. 



Inserting ^ into (O and eliminating time leads to a 
set of equations for the amplitudes Bkn- These equa- 
tions yield a characteristic amplitude of the interaction 
between the open channel Bon and the closed ones Bfe^o.n 
inside the DB core, which will be denoted by Vac-dc- If 
we approximate -B/t^o.n — Oi coefficients in the remaining 
equations for _Bon describe the time-averaged or dc part 
of the scattering potential. 



DETERMINING THE LOCATION OF 
RESONANCES 

As was shown above, a breather acts as a time periodic 
scattering potential for propagating waves. The system 
© with time periodic coefficients with period Tt can be 
considered as a Floquet problem. This system possesses 
stable solutions, which satisfy the following condition 



(11) 



These are so-called Floquet or Bloch eigenstates. Here 
is a Floquet multiplier, which can be represented as 
6 = ujTb with some frequency to. Note that, strictly 
speaking, such frequencies are defined only modulo flh- 
Most of all Floquet states are extended ones with fre- 
quencies to belonging to the spectrum wi.2,3('?)- But due 
to the spatial localization of a DB, some Floquet states 
could be also localized around a DB with frequencies uj 
corresponding to some localized modes lul. 

In order to predict the location of a Fano resonance 
we use the results of [ij ^7 19]- The resonance posi- 
tion is determined by the bound states of the dc part of 
the scattering potential if the width of the resonance is 
small compared to the continuum band width. That con- 
dition is equivalent to requesting that the squared width 
1//3| of w|((7) is larger than V^^^c-dc- introducing the 
parameter 



A, 



(12) 



it follows Xf 1- In this case the bound states of the 
dc part of the scattering potential ujf'^ can be considered 
as additional discrete levels, which are weakly but reso- 
nantly coupled to the open channel via the frequency of 
the DB ilft. In other words, the Fano resonance takes 
place when 



dc 



(13) 



is satisfied for some qp- In order to predict the location 
of the resonance we have to know the frequencies of the 
bound states of the dc part of the scattering potential. 
We obtain these values numerically by diagonalizing the 
corresponding matrix. 

Let us estimate Vac-dc for our case. In the breather 



core the cos(ribi) = 0.5(e 



-iQht 



) terms in (O yield 
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a coupling between channels k and A: ± 1 of the order 0.5. 
Thus Vac-dc ~ 0.5. Consequently we obtain Xp = 0.5/3|. 
Thus for < 1 we can predict the position of the Fano 
resonance as described above. Our strategy of searching 
for Fano resonances is to compute the frequencies of lo- 
calized modes of the dc part of the scattering potential 
and to satisfy (|13|l . 



COMPUTATIONAL RESULTS 

For our computations of the transmission coefficient 
we studied two types of breathers with left-right (A-type) 
and up-down (B-type) symmetries with one vertical re- 
sistive junction (see Fig. |2J|. 



two of them, which satisfy the condition H13|l : w^'i ^ = 

0.62, wf2,A = 0.73, ljII b = 0.93 and wfa.s = 0-95. ' 

Thus we expect to observe two Fano resonances for the 
A-type breather at frequencies Uq^ = 3.72 and Uq^ = 
3.83 according to H13|) and two for the B-type breather 
at frequencies ujqp = 4.03 and ujqp = 4.05. Indeed, the 
computation of the transmission coefficient shows that 
there are two Fano resonances for the A-type and B-type 
breathers (Fig. ^J. The frequencies of the numerically 
observed resonances are iVq^ — 3.68 and ujqp = 3.81 for 
the A-type breather and ujqp — 3.93 and ujqp = 4.07 for 
the B-type breather (see Fig. ^ . 



The Floquet approach 

We use the formalism described above, and calcu- 
late the transmission coefficient in the Hamiltonian limit 
was by implementing the numerical scheme described in 
Ref. . According to that scheme, we use proper bound- 
ary conditions in the form of a time-periodic drive with 
a given frequency from the spectrum l|21|) . We simulate 
a system with = 20 cells. On the right hand side of 
the system we use an additional restriction which implies 
that there is exactly one transmitted plane wave, propa- 
gating to the right. On the left hand side we do not apply 
any additional restrictions. The local time-periodic drive 
generates a mixture of waves with different amplitudes 
propagating to the left and to the right. In general, we 
excite a superposition of all possible polarization vectors. 
However, we drive our system with a frequency from one 
dispersion curve. It means, that all waves which cor- 
respond to other dispersion curves decay in space and 
the desired scattering setup is realized at some distance 
from the left boundary. In our case the two additional 
branches uji^2{<l) are dispersionless and their waves do 
not propagate through the system. Therefore, we can 
already choose the first cell away from the boundary in 
order to successfully apply the numerical scheme from 
Ref. 13]. If we take into account a nonzero damping, 
then the seond band uj2 (q) becomes weakly dispersive as 
well, and evanescent modes at these frequencies pene- 
trate slightly into the ladder. We then simply move our 
reference sites for the scattering setup further away from 
the edges of the ladder, where again a single frequency 
excitation is found. 

For Pl = 0.5 and 77 = 0.35 the spectrum (|21|) is located 
between uj3{0) = 3.53 to tt'3(7r) = 4.52. We have found 
by numerical diagonalization that for f2f, = 3.1 all local 
modes of the dc part of the scattering potential of A- and 
B-type breathers with one vertical resistive junction are 
located below the plasma frequency uji and there are only 
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FIG. 4: Transmission coefficient versus 0^3(5) with Fano res- 
onances for waves scattered by A-type (dashed) and B-type 
(solid) DBs with one vertical resistive junction for Qt = 3.1, 
/3l = 0.5 and r; = 0.35. 

The number of localized modes of the dc part of the 
scattering potential does not change when the damping 
and bias become nonzero. But in this case there are 
two dispersion curves with nonzero dispersion 0^2,3 (9)- 
One of them is located above the plasma frequency 
(jJi(q)l < \uj3{q)\ (as before) and the second one below it 
^2(9)1 < |'^i(9)|- Since the latter one possesses a rather 
weak dispersion, for small damping a ~ 0.01 there are 
still localized modes, which satisfy the condition (|13ll . 

We also incorporated the finite damping into the calcu- 
lation of the transmission coefficient by using the trans- 
formation The resulting scattering problem can be 
again analyzed along the lines of the Floquet approach 
from above. Note that this modified scheme implies the 
excitation of a scattering wave setup which decays in time 
in a spatial homogeneous way. The numerical results of 
this scheme for a = 0.01 coincide with the transmission 
shown in Fig.^ The Fano resonance positions are practi- 
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cally the same, because the frequencies of locahzed modes 
almost did not change. 

Direct numerical simulations 

In order to test whether the Fano resonances com- 
puted along the lines of the Floquet approach (see Fig.^ 
also persist in an experimental situation, we performed 
direct numerical simulations of the scattering of linear 
waves by a DB in Josephson ladders which emulate a 
real experimental situation. The equations of motion 
were numerically integrated using a standard fifth or- 
der Runge Kutta scheme |22| with precision monitoring. 
To study the scattering of linear waves by a DB in a 
situation comparable to an experiment, a special scat- 
tering experiment scheme was set up. Linear waves are 
generated in a JJL with open ends by locally applying 
a time-periodic current 71 (t) at the vertical junction 1: 
71 (t) = jac cos(wi). The local current acts as a local para- 
metric drive. It excites a tail of junctions that oscillate 
with frequency lu. This tail extends into the ladder and 
decays exponentially in space. The situation is compa- 
rable to a DB with fre quen cy oj located at the system 
boundary. In Refs. [isl l2nl|. the spatial extent of oscil- 
latory tails of DBs was computed analytically. If the 
frequency uj lies inside the dispersive L03 band, the tail 
decays as ipn oc exp(— an). Outside the band, the oscil- 
lations practically do not penetrate into the system. To 
monitor the linear wave propagation in the system, we 
compute the time-averaged oscillation power of vertical 

junction n, Pac.n = iVn)- 

As a scatterer, a DB of either type A or B is launched 
in the center of the system. The DB frequency is then 
tuned by a spatially uniform DC bias 7. 

We simulate an open-ended JJL with N — 100 cells, 
damping a = 0.05, discreteness (3l = 0.5, and anisotropy 
77 — 0.35. The vertical junction at the left edge n — 1 
is driven by a time-dependent current with 7ac = 0.05. 
The driving frequency is swept within 3 < w < 5. At 
each frequency step, the system is integrated over a time 
period of 500 units to allow for relaxation. The dynami- 
cal values are then monitored over the following 300 time 
units. During that time the oscillation power and phase 
velocities are averaged and recorded. 

Fig. Ufa) displays the obtained average oscillation 
power spectrum in the vertical junction n = 101 (at 
the right edge of the ladder which is opposite to the ac- 
driven junction). In the empty case without a scattering 
DB state inside the ladder, the oscillation power shows 
a broad peak when the excitation frequency uj is inside 
the UJ3 band. Outside the band, the oscillations quickly 
decay to a level below 10"^*^. As a test, the frequency- 
dependent spatial decay was compared to the analytical 
predictions and perfect agreement was found. When a 
scattering DB is inserted into the system at the central 



site n = 50, the propagation of linear waves is modified, 
and dips may appear in the power spectrum (dashed line 
in Fig. El 

We obtain the transmission coefficient r by relating 
the oscillation power at the boundary vertical junction 
(n — 101) with and without a DB inside the system, as 

fac,ioi(with DB) 
~ Pac.ioi (without DB)' ^ ' 

Fig. |SJb) shows the obtained transmission for type A 
scatterer states. Outside the 0J3 band, linear waves do not 
propagate, hence any obtained transmission coefficient is 
meaningless. Inside the band, we show the transmission 
through DBs of three different frequencies. In the first 
case of small DB frequency fib = 2.111 ^ the transmis- 
sion curve has a broad maximum. The DB is practically 
transparent in the center of the band. If we increase 
the frequency of the scattering DB, dips appear in the 
transmission at frequencies Wdip. The values of these 
frequencies change upon variation of the scatterer DB 
frequency flh- 

We systematically varied the scattering DB frequency 
and computed various transmission curves. In Fig. [SJ 
we show the dependence of the positions of the dips in 
the spectrum w^ip as a function of the scattering DB 
frequency fl^ for type A and type B states. The obtained 
dependence is linear and follows 

t^dip = Ao + Aif^fa. (15) 

For the type A scattering DB, two dips were observed. 
For the lower one we find Ag = 0.601 and A-y = 1.015 
from a least squares fit, while the upper resonance satis- 
fies Aq = 0.680 and A\ = 1.009. In the case of a B type 
DB, we observe only one resonance, following Aq = 0.995 
and Ai = 0.966. We interpret Aq as the frequency of a 
dc local mode, while Ai is practically equal to unity, and 
find perfect agreement with Eq. fl^ . We thus conclude 
that the dips in the transmission spectrum are in fact 
Fano resonances. 

The persistence of the Fano resonance in the direct 
simulations of the finite size Josephson ladder including 
damping gives a strong indication for a possible experi- 
mental observation of the phenomenon. Experimentally, 
the time-dependent local biasing of one ladder end re- 
quires microwaves. However, instead of using an external 
source, a resonant DB state, as studied experimentally in 
Ref. |21j, could be used to create hnear waves. A scatter- 
ing setup should therefore consist of a resonant DB state 
used as an emitter (creating linear waves) , while another 
(non-resonant) DB within a suitable distance would act 
as a scatterer. For the detection of radiation we suggest 
the use of another set of resistive Josephson junctions 
inside the ladder, situated at the far end. These junc- 
tions may be seen again as a DB state. Their property of 
quasi-particle detection can be then used by biasing them 
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frequency co frequency co 

FIG. 5: Direct numerical simulation of the linear wave propagation in a JJL with A'^ = 100 cells, a = 0.05, /3l = 0.5, rj — 0.35 
with a boundary ac bias 71 — 7ac cos(tjt). (a) Oscillation power of vertical junction n = 101 versus excitation frequency to for 
an empty system (solid grey line), for an A DB at site n = 50 with frequency fit = 2.777 (solid line), and fit = 3.284 (dashed 
line), (b) Transmission coefficient t for an A DB of frequency Qi, = 2.777 (solid line), fi;, = 3.284 (dashed line), and frequency 
n = 3.082 (dotted line). 
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FIG. 6: Dip frequency ujdip in the transmission from Fig.|3^b) 
versus DB frequency for a type A DB (two dips, marked by 
circfes and squares) and a type B DB (one dip, diamonds). 
Lines indicate linear fits of the dip frequency to Eq. II15II . 



close to the superconducting gap voltage '23l|. Fig.[7|gives 
a sketch of the experimental setup of locally biased DB 
states used as linear wave source, scatterer, and detec- 
tor. Promising experimental studies are under way (see 
Rcf. 24] for details). 



CONCLUSION 

We investigate theoretically the existence of Fano res- 
onances in wave scattering by DBs in JJLs. Due to the 



weak coupling between open and closed channels of the 
scattering potential, generated by a DB, perfect reflection 
occurs, according to the theory in the Hamiltonian limit, 
when multiples of frequency of DBs match the difference 
between the frequency of a dc local mode and some fre- 
quency in the spectrum. Numerical calculations of the 
transmission coefficient based on the Floquet approach 
confirm this analytical prediction. Moreover, Fano reso- 
nances survive even the presence of dissipation and exter- 
nal dc bias. Direct numerical simulations of the damped 
and biased system demonstrate Fano resonances as well. 
This is a strong indication that Fano resonances can be 
observed experimentally in Josephson junction ladders. 

An interesting question is why the direct numerical 
simulations do not show a perfect reflection (cf. Fig.[SJ), 
at variance with the computation of the transmission co- 
efficient based on the Floquet approach (cf. Fig.^. We 
checked that this is neither an artifact of the numerics, 
nor due to finite size effects. The two methods (direct 
simulations and Floquet based approach) differ in the 
way boundary conditions are defined. The only explana- 
tion we thus arrive at, is that although the damping is 
weak, it leads to a different damping of the propagating 
waves in different channels. Waves propagating locally 
around the breather at different frequencies (i.e. in dif- 
ferent channels) are damped at a different rate, which de- 
pends on the frequency of oscillations. Correspondingly 
the phase and amplitude relations between the partial 
waves in different channels are changed, leading only to 
a partial destructive interference . Alternative expla- 
nations involve the broadening of a resonance frequency 
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FIG. 7; Schematic experimental setup for measuring the plasmon scattering by a DB, showing source, scattering, and detector 
DB. The individual DBs are controlled using local bias currents 71-73. Technical details are found in Ref. |24j]. 



line due to damping. The damping causes the trans- 
mission to stay nonzero in the resonance, and changes 
(broadens) the line width of the resonance. This also 
leads to the impossibility of resolving two closely nearby 
lying resonances, exactly as we found in our numerical 
studies for the type B breather. 

Fano resonances can be considered as a benchmark 
of dynamical localized excitations (DBs) in their action 
on propagating waves in the system. A similar proposal 
for the observation of Fano resonances in light-light 
scattering in nonlinear optical media has been reported 
in Ref.Hg. 
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APPENDIX A: DISPERSION OF LINEAR WAVES 

Substituting the expressions ^ into system to- 
gether with the transformation (PJ and linearizing with 
respect to the small excitations, we obtain 



a 
T 



1 

Pi 



+ = -^i^<+^n-^n) (16) 

4 tiPl 



(1 



1 



4'^" 11(3 L 
By using the plane wave ansatz 

= e'(«"-"'*)i*, , 

one obtains that there are three bands, 
dispersionless 

9 

u;l{q) = 1 - — . 



(17) 

One of them is 
(18) 



For a = this is the plasma frequency and it is character- 
ized by in-phase excitations of upper and lower horizontal 
junctions = and all vertical junctions being at rest 



A^ — Q. Due to the absence of dispersion this mode can 
be excited in each cell with arbitrary amplitudes, since 
it does not propagate along the ladder. 

The two other bands are located above and below it 



F = 



G 



111 / 1 



2 
(1 



2 



72 + -^(1 -cosg) - ^ 



7 

' 4 ' 



(19) 



2 ,^ , a^, a^, 
--(l-(l + -)cos,--) 



The dispersion of the second band LU2{q) is weak. The 
polarization vectors of these two branches can be written 
in a compact form 



A'' = -A'' , A'' = 



1 



A" (20) 



In the Hamiltonian limit a = and 7 = the band 
102(1) — 1 becomes dispersionless as well. Its polariza- 
tion vectors are simplified to out-of-phase excitations of 
upper and lower horizontal junctions A'^ — —Ah and 
A'' = i(l - e'9)A". The third band ^3(9) keeps a fi- 
nite dispersion 



^3(9) 



1 



2 2 



cos q) 



(21) 



with corresponding in-phase excitations of upper and 
lower horizontal junctions and the following relation be- 
tween vertical and horizontal ones 



A"" = A" 



A" = 77(e-'' - l)Ah 



(22) 
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